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ABSTRACT 

A spectroscopic analysis of SDSS J160043. 6+074802. 9, a binary system containing 
a pulsating subdwarf-0 (sdO) star with a late-type companion, yields Toff — 70 000 
^ ■ ± 5000 K and log g = 5.25 ± 0.30, together with a most likely type of K3 V for 

\ the secondary sta r. We compare our results with atmospheric parameters derived 

Cf^ , by I Fontaine et al.l ()2008D and in the context of existing evolution models for sdO 

0^ ' stars. New and more extensive photometry is also presented which recovers most, 

\ but not all, frequencies found in an earlier paper. It therefore seems probable that 

. some pulsation modes have variable amplitudes. A non-adiabatic pulsation analysis of 

' uniform mctallicity sdO models show those having logg > 5.3 to be more likely to be 

^\ , unstable and capable of driving pulsation in the observed frequency range. 
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1 INTRODUCTION Pickering lines of He II. No He I lines (such as the typ- 
ically strong features near 4471 and 4026A) are detected, 

In the course of a search for new AM CVn stars, j-uUng out types sdB and sdOB. The relatively narrow ab- 

the Sloan Digital Sky Survey (SDSS) object, SDSS sorption hnes also seem to rule out the possibility of the 

J1600 43.6+074802.9, wa s discovered to be a very rapid pul- object being a white dwarf star. In the red region, the Na D 

sator IWoudt et alj|200^, hereafter Paper I). Subsequent ob- Unes (5889 and 5896A) are present and - given the presence 

servations, reported in the same paper, shows the variability of He II lines - indicate that the object is a binary. This 

to be complex, consisting of a strong variation (amplitude ^as confirmed by changes in radial velocities determined 

-0.04 mag) with a period of 119.33 s, its first harmonic from two spectrograms obtained a few weeks apart, though 

(0.005 mag; 59.66 s) and at least another eight having peri- the data were insufficient to even suggest an orbital period, 

ods between 62 and 118 s and amplitudes in the range 0.003 The spectru m is apparently that of a classical sdO star (see, 

to 0.007 mag. for example. iMoehler et al]ll990l ) combined with a late-type 

Preliminary spectroscopy in Paper I showed the star Main-Sequence star, 

to be a hot (sdO) subdwarf with a cooler Main-Sequence For convenience, we abbreviate the object name, SDSS 

companion. The object was known to be very hot from the J160043.6+074802.9, to J16007+0748. The 2000 equatorial 

SDSS colours («' — g' = —0.201, g' ~ r' = —0.205) and the co-ordinates are implicit in the full SDSS name; the equiv- 

spectrograms confirm this, showing He II 4686A and the alent galactic co-ordinates are I = 18.65°; 6 = +41.27°. The 

Pickering series lines of He II at 5411, 4542 and 4200A to Sloan Digital Sky Survey gives photometry: u — 17.21, g' = 

be clearly present in absorption. The Balmer series lines H/3 17.41, r' — 17.61, i' — 17.66, z' = 17.77; for more details and 

to He are also present, but these are probably blended with background, see Paper I. 
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In this paper we present a spectroscopic determination 
of effective temperature (Teff), surface gravity (logg) and 
helium abundance {N{He)/N{H)) using a flux calibrated 
SDSS spectrum. We also present new and more extensive 
photometry for Jf6007+0748, which recovers most, but not 
all frequencies found in Paper I; we conclude that some 
of the pulsation modes have variable amplitudes. Several 
uniform metallicity sdO star models, having Toff and log g 
within the 68 per cent confidence limits derived for these 
quantities, were tested for non-adiabatic pulsation stability. 
However, pulsati ons are not truly driven in our models, in 
contrast with Fo ntaine et al.l (|2008l ) non-uniform metallicity 
sdO models, which achieve real destabilisation. 



2 SPECTROSCOPIC ANALYSIS 

In this section, we describe an analysis of a SDSS spectro- 
gram of the system taken from the sixth data release of th e 
Sloan Digital Sky Survey (lAdelmann-McCarthv et al.ll2008l ). 

A late-type companion (Paper I) with a luminosity 
high enough to produce observable Na I D lines, and in- 
sufficient spectra with which to d etermine an orbi t al pe - 
riod, meant that the procedure of I Simon fc SturmI (|1994| ) 
could not be used to disentangle contributions made by 
the two stars at each wavelength. Because a flux cali- 
brated spectrum of J 16007-1-0748 had been obtained on 
2006 April 28th in excellent weat h er, a s part of the 
SDSS survey, the lO'Donoghue et al.l l|l997l ) technique for 
analysing energy distributions of composite sdB stars was 
developed and applied. Specifically, as Paper I identifies 
J 16007-1-0 748 as a sdO and late- t ype Main Sequence bi- 
nary, the lO'Donoghue et al.1 (119971 ) method was extended 
to simultaneously extract all parameters of interest from 
the SDSS energy distribution; these are the interstellar 
reddening, Tefj, angular radii (6) and radial velocities (v) 
of both stars, as well as the logg and helium abundance 
{N{He)/N{H)) of the sdO star. 

The program tlusty (|Hubenv fc Lanz 1995) was used 
to compute a grid of model stellar atmospheres for the sdO 
in Non-Local Thermodynamic Equilibrium (non-LTE) , com- 
prising hydrogen and helium only; this was justified only in 
so far as lines attributable to other elements, and originat- 
ing in the sdO star spectrum, are not identified in Paper I 
spectra. Standard model atoms for H I, He I and He II were 
adopted with the number of non-LTE levels being 9, 14 and 
14 respectively. Level populations for all non-LTE levels were 
calculate d with occupation probabil i ties t aken into account 
following iHubenv. Hummer fc Land (|l994l ): this also applies 
to the calculation of synthetic spectra described below. 

With computed model atmospheres, non-LTE synthetic 
spectra (3500 A ^ 68OOA and having hydrogen and he- 
lium lines only) were calculat ed using th e program syn- 
SPBC by Hubeny & Lanz (see IZborill Il996l . and references 
th erein) with revisions to He I line broadening as described 
bv lMortimore fc Lvnas-Gravl l|2006l ). Non-LTE synthetic hy- 
drogen and helium line spectra were produced for the sdO 
star for the wavelength range 3500 A ^ 68OOA. Balmer 
line b roadening was based on the Vidal. Cooper fc SmithI 
l|l970h static ion approximation as ISchonind l| 19941 ) shows 
ion dynamics mainly affect line centres and are masked in 
stellar atmospheres. He II line broadening was calculated 



using tables by ISchoning fc Butler! l|l989l ). For comparison 
with observation, theoretical spectra were convolved with 
the appropriate wavelength dependent instrumental profile; 
for reasons discussed below, the projected stellar rotation 
velocity (v sin i) was ta ken to be zero. Inte rstellar reddenin g 
was modelled foUowing lSeatonI (| 19791 ) and lHowarthI l|l983l ). 

An initial manual comparison of observed and theoret- 
ical spectra was based on the grid 55 000 ^ T^s ^ 70 000 K 
(AToff = 5000), 5.0 ^ log 3 < 6.5 (A logs = 0.5) and 
0.1 ^ N(He)/N{H) <^ 0.5 {A{N{He)/N{Hy) = 0.2); from 
this it became clear the sdO star had a T^s near or above 
the upper limit of the initial grid. After some trials, the fi- 
nal search was based on linear interpolation in the parameter 
space cube having {T^s ,\og g , N [H e) / N [H)) coordinate ex- 
tremes specified by (66 000 K,4.9,0.2) and (76 000 K,6.5,0.5). 

For the late-type companion, synthetic spectra by 
iMartins et al.l (2005) were used, in preference to atlases de- 
rived from observed spectra, because there is then a well 
defined T^s scale on which to base interpolation. Paper I 
results suggest the companion is a late F or early G-type 
Main Sequence star. Accordingly, synthetic spectra for the 
late-type companion were generated by piecewise linear in- 
terpolation in the range 3000 ^ Tcff ^ 9500 K; loggf = 4.5 
and solar metallicity being adopted in the first instance. 
Again the projected stellar rotation velocity was assumed 
to be zero and synthetic spectra were convolved with the 
appropriate wavelength dependent instrumental profile. 

Trial radial velocities between —200 and +200 km/s 
were selected for both components of the binary. The max- 
imum fiux density among all sdO star model spectra was 
identified, and divided by the maximum fiux density in the 
observed energy distribution; this ratio when multiplied by 
10.0 corresponded to the maximum trial angular radius for 
the sdO star. A maximum trial angular radius for the late- 
type companion was similarly selected, except that the al- 
most arbitrary scaling factor chosen in this case was 50.0 to 
refiect the smaller contribution to the observed energy dis- 
tribution and the longer wavelength of the peak flux density 
in its energy distribution. Zero was selected as the minimum 
trial angular radius for both stars and trial interstellar red- 
dening values were in the range 0.0 ^ Eb~v ^ 0.4. Because 
interstellar reddening may be signiflcant, and associated in- 
terstellar Ca II and Na I D lines present in the data, the 
Ca II (K) line and central pixels of both Na I D lines were 
ignored during the fitting process. 

Shifting and scaling appropriately chosen theoretical 
spectra by correct radial velocities and squared angular radii 
respectively, then adding them together and adjusting for 
interstellar reddening, should reproduce the observed en- 
ergy distribution; this fitting process was carried out us- 
ing a revision (Vers ion 1.2) of the genetic algorithm which 
ICharb onnea ul (|l995l ) codes as pikaia. Briefiy, nine unknown 
parameters were regarded as the "genetic signature" of an 
individual whose fitness was determined by how well the 
observed energy distribution was reproduced; this was char- 
acterised by the inverse chi-squared (x~^)- The principle of 
natural selection means fitter individuals are more likely to 
breed and produce fitter offspring. Continuing evolution of 
a population of individuals over many generations eventu- 
ally locates a global minimum which represents the best 
obtainable fit. 

On the basis of constraints discussed above, a random 
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Figure 1. Observed SDSS energy distribution and estimated con- 
tributions from the binary components. 



number generator was used to create populations of indi- 
viduals. Initially populations of 100 individuals, selected 
with 100 different random number generator seeds, would 
be evolved for a 100 generations; those that evolved most 
rapidly would be selected for further runs, evolving them 
over larger numbers of generations. The global minimum 
was identified as the almost unique reached after evolv- 
ing three distinct populations for 30 000 generations. Further 
tests, varying metallicity and surface gravity of the late-type 
companion gave less satisfactory fits. Fig. 1 shows the ob- 
served SDSS energy distribution and deduced contributions 

to it from both binary components. 

Standard deviations were estimated following lAvnil 

l|l976t l. having regard for the simultaneous estimation of nine 
parameters. Hypersurfaces of equal in nine-dimensional 
parameter space were constructed from all trial solutions 
generated during an evolution of 30 000 generations. For nine 
degrees of freedom, one standard deviation corresponds to 
A(x^) = 10.5. Parameter values and standard deviations 
are presented in Table 1, where it should be noted that he- 
liocentric corrections have not been applied to the derived 
radial velocities. 

In Fig. 2, a thin line shows the observed SDSS energy 
distribution compared with the sum of the sdO and late- 
type Main Sequence energy distributions (shown in Fig. 1) 
plotted as a thicker line; the "spike" near 5600A has been 
removed. Note that the strong Ca II K-line near 3933 A is 
not fitted and He is predicted to be weaker than actually 
observed. There is clearly an interstellar contribution but 
also note that Ca II (K) appears anomalously strong in the 
SDSS spectrum when compared with the Paper I (fig. 6) 
spectrum. 

Otherwise the fit shown in Fig. 2 appears to be satisfac- 
tory. In particular, features near 5100 A attributable to the 
late- type companion are reproduced and justify Tcff (K3 V) 
in Table 1 rather than the higher value (corresponding to 
a late F or early G Main Sequence star) which Paper I 
suggests. Moreover there appears to be little prospect of 
obtaining a superior fit to the SDSS energy distribution us- 
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Figure 2. Fits to energy distribution of J16007-f 0748. Central 
wavelengths, corresponding to a displacement of O.oA are entered 
in each panel. 



Table 1. Derived Parameter Values and Standard Deviations 



rcfi(sdO) = 70 000 ± 5000 K 
logg(sdO) = 5.25 ± 0.30 
N{He)/N{H){sdO) = 0.31 ± 0.11 
e(sdO) = 6.27 ±0.25 x IQ-^^ radians 
i;(sdO) = -66 ± 26 km/s 
Teff(K3V) = 4690 ±75 K 
6I(K3V) = 5.26 ± 0.15 x IQ-^^ radians 
i;(K3V) = -67 ± 32 km/s 
Eb-v = 0.173 ±0.013 



ing theoretical profiles additionally broadened as a conse- 
quence of stellar rotation; for this reason, the assumption of 
II sin i — 0.0 for both stars was appropriate. 

For 2326 degrees of freedom (the number of wave- 
length points used less the number of free parameters), 
^/l^x^) = 67.88 was obtained; statistically y^2x^) should 
approximate to a Gaussian distribution having a mean of 
68.21 and unit variance. Counts in adjacents spectrum pix- 
els are correlated and so a good fit may not be verified using 
the x^-test. With the Runs Test (|Wald fc Wolfowit j [l940h 
the expected number (< r >) of alternating sequences of 
positive and negative fit residuals was estimated. For the 
Fig. 2 fit, with 2335 valid data points, < r >= 1166.5 ±24.1 
(s.d.); the observed number of runs was r — 914 which is 
significantly smaller than < r >, as may be expected (in 
part at least) having used H/He synthetic spectra to rep- 
resent the sdO star which would be expected to have weak 
(unobserved) metal lines. 
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3 INTERSTELLAR REDENNING 

In order to check the Section 2 solution, and determina- 
tion of Eb-v in particular, populations of 100 individuals 
were evolved for 10 000 generations having Eb-v fixed at 
0.0, 0.044 and 0.173. Eb-v = 0.044 was chosen because 
it is predicted by the adopted galactic extinction model as 
explained below. All other parameters were free, with the 
exception of T^s , which was fixed at 70 000 K for reasons 
discussed below, with values determined by attempting to 
fit the observed energy distribution as already described. For 
Eb-v = 0.0, 0.044 and 0.173 the values of a/(2x^) obtained 
were 74.87, 71.81 and 67.97 respectively; these are 6.66, 3.60 
and -0.23 standard deviations away from the expected 68.21. 
As already noted, individual fluxes in adjacent pixels of the 
SDSS spectrum are not independent and so Eb-v = 0.0 or 
0.044 cannot be ruled out with a high degree of confidence; 
although it was clear a satisfactory fit could not be obtained 
with these interstellar extinctions. 

The observed energy distribution is strongly dependent 
on Eb-v and, for a hot star such as J16007+0748, only 
weakly dependent on TcS- reff(sdO) and log(j(sdO) are ob- 
tained from simultaneously fitting available Balmer lines. 
Forcing Eb-v = 0.0 or 0.044, with rcfi(sdO) unconstrained, 
inevitably drove solutions towards regions of parameter 
space where refi(sdO) was unphysically low, given the ob- 
served spectrum. Of particular interest is using Eb-v = 0.0 
and reff(sdO) = 70 000 K, as it all ows a direct compar- 
ison with the iFontaine et al.l (|2008l ) result. We obtained 
\ogg ~ 5.37 ± 0.30, which is still off by 0.15 dex of their 
lowest \ogg determination. 

Interpretation of the derived Eb-v requires a knowl- 
edge of the wavelength-dependent light loss (if any) at the 
spectrograph slit while the spectrum of J 16007-1-0748 and 
any flux calibration standards used by the SDSS project, 
were being obtained. Bohlin fc Gilliland's (2004) absolute 
flux distribution was used with lSmith et aL (120021') u ' n'r'i'z ' 
filter response curves, and equation 3 in lSmith et ahl i|2002l ) . 
to obtain synthetic magnitudes for the SDSS standard 
BD -|-17°4708; these agreed with photometric results within 
published error limits. Applying the same approach to 
the SDSS flux cahbrated spectrum of J16007-f0748 gave 
((?' — r') = —0.176 for comparison with the photometric 

value of -0.205. 

Applying the transformation given by ISmith et al.l 

l|2002h (their table 7), gave the corresponding colour cor- 
rection, due to wavelength-dependent light loss at the spec- 
trograph slit, of A(B - V) = 0.030. Reddening in the 
direction of J16007-I-0748 was therefore estimated to be 
Eb-v = 0.143. Of course, the assumption made here was 
that J16007+0748 only varies as a consequence of the pul- 
sation Paper I reports. 

For a high galactic latitude object hke J16007-I-0748, 
the axisymmetric mo del (which does n o t incl ude a spiral 
arm contribution) by lAmores fc Lepind (|2005h was appro- 
priate for estimating the galactic contribution to interstellar 
reddening along its line of sight; this gave Eb-v = 0.044 
for all distances greater than 500 pc. The contribution to 
Eb-v unaccounted for was then 0.099. Anomalous redden- 
ing in the direction of J16007-I-0748 is one possibility, as is 
the presence of dust and gas in the vicinity of this object 
which could have come from the sdO progenitor; these and 



O 



5.0 



5.5 



6.0 



6.5 




Figure 3. Predicted positions of J16007+0748 in a Teff-log g dia- 
gram. The large upper rectangle shows the position and its uncer- 
tainty as determined in the present paper. The lo wer and smaller 
rectan gle is the corresponding position which iFontaine et al.l 
1I2OO8I) determine. See text for further details. 



other possibilities need a further investigation, which is be- 
yond the scope of the present paper. 



4 EVOLUTION 

While work reported in this paper was nearing completion, 
IFontaine et al.l (|2008l ) published their spectroscopic analysis 
of J16007-I-0748 and explain Paper I pulsation in terms of 
radiative levitation. The log(;(sdO) Fontaine et al. obtain is 
5.93 ± 0.11 which disagrees with our Table 1 result, though 
the corresponding T^s values are consistent within error lim- 
its. The weakness of our analysis was the dependence on one 
SDSS spectrum, an d the SDSS flux calibr ation in particular. 
The strength of the lFontaine et all (|2008l ) work lies in those 
authors having obtained a high signal-to-noise spectrum, al- 
though of lower resolution than the SDSS spectrum used in 
the present work, on which to base their analysis. 

A sequence of Balmer lines (Ha - H8) provides a Te ff 
and log g determination as is well known (jSaffer et al .11 1993 ) ; 
this fact allowed authors of the present paper, using the 
observed energy distribution, to simultaneously determine 
Eb-v (qualified as above) which Fontaine et al. implicitly 
adopt as zero. On the other hand, although the Fontaine 
et al. subtraction of a late-type template spectrum is not 
quite perfect (notice the residual Na I D line in the sdO 
star spectrum in their fig. 1), they do identify the late- 
type companion as a GO V star which is consistent with 
Paper I. Giv en the sim i lar T eff(sdO) determinations ob- 
tained bv F ontaine et al] (I2OO8I ) and ourselves, the account- 
ing for interstellar reddening in our study required a redder 
model of the observed energy distribution; this was achieved 
through a redder late-type companion, which might qualita- 
tively account for the difference between the two analyses. 

Determinations of Tcff and logg for J16 007+0748 are 
compa red with those for 46 sdO stars by IStroeer et al.l 
(|2007l ) in Fig. 3. Helium-deficient sdO stars are plotted 
as and helium- e nrich ed sdO stars as"*". Uncertain- 

ties in IStroeer et all (|2007l ) measurements are represented 
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by a cross plotted in the top left-hand corner. Also shown 
in Fig. 3 as solid lines are three Post Asymptotic Giant 
Branch (AGB) tracks bv lBlockeil l|l995al lbh: different tracks 
are distinguished by different colours (blue - hydrogen burn- 
ing 0.605 Mq, green - helium burning 0.625 Mq and red - 
helium burning 0.524 Mq). Dashed li nes are used to plot 
helium white dwarf evolution tracks bv lDriebe^ et all l|l998l . 
I1999I ) and the position of the Zero-Age Helium Main Se- 
quence shown with a dotted line. Each evolution track plot- 
ted in Fig. 3 is lab el led wi th the corresponding stellar mass. 

IStroeer et all (|2007l ) discuss evolutionary scenar- 
ios for both helium-deficient and helium-enriched sdO 
stars. J16007-I - 0748 is found to be helium-enriched by 
iFontaine et al.l ()2008l ) and confirmed, in principle at least, 
by results in Table 1. Only Stroeer et al. scenarios concerned 
with helium-enriched sdO evolution were therefore consid- 
ered^ 

iDriebe et al.l l| 19981 ) evolve 1 Mq stars from pre-Main 
Sequence to the Red Giant Branch (RGB); their envelopes 
are almost completely stripped by Roche-lobe overflow in a 
close binary system before core helium burning starts. The 
remaining hydrogen-shell turns out to be rather thick and 
it is not obvious how a helium-enriched atmosphere is ob- 
tained. If the r esidual mass li es in the range 0.21 Mq ^ 
M ^ 0.3 Mq, iDriebe et all (jl998h show that hydrogen- 
shell flashes occur which may lead to a dredge-up of he- 
lium, but predicted positions of J16007-I-0748 in Fig. 3 are 
inconsistent with Driebe et al. tracks having these masses. 
If J16007-I-0748 is a post-RGB object, it has to have a mass 
~ 0.4 Mq and it is not then clear how helium-enrichment 
could have occurred. Continued stripping of the hydrogen 
envelope by interaction with the late-type companion is a 
possibility though this would have speeded up evolution, 
because most of the luminosity comes from the hydrogen 
burning shell, making detectability less likely. 

Positions of J16007-)-0748 ir i Fig. 3 are also c l ose to 
post- AGB tracks (not shown) bv ISchonberneJ (|l979l . 1198^ ) 
and, like the heliu m-enriched sdO HE 1430-0815 which 
IStroeer et all |20o3) consider, it could be a post-AGB star, 
but the timescale for evolution from the AGB to the pre- 
white dwarf stage l|Schonberneijll97^ 1 is only ~ 3 x 10* years 
and again this reduces the probability of de tection. Unfor- 
tunate ly the pulsation mass of 0.48 Mq by iFontaine et al.l 
l|2008t ). which these authors choose to be consistent with 
typical sdB masses and therefore needs improvement as 
they acknowledge, did not enable us to distinguish between 
the post-RGB (M ~ 0.4 Mq) and post-AGB scenarios 
(M ~0.55Mq). 

iBlockeij fmEaV) develops the earlier work by 
ISchonberneil (|1979, . ,19831 ) by considering evolution of low 
and intermediate mass stars with mass-loss taken into ac- 
count with different descriptions for the RGB, AGB and 
post-AGB phases. Three of Blocker's post-AGB tracks pass 
close the error rectangles in Fig. 3. The hydrogen burn- 
ing 0.605 Mq track (shown in blue), for which the pro- 
genitor mass was 3 Mq, passes near the top left-hand cor- 
ner of the error rectangle obtained in the present paper, 
which suggested the sdO in J 160074-0748 would have to have 
M ^ 0.605Mq if it is a hydrogen burning post-AGB star. 
The post-AGB age would be only ~ lO** years and detecting 
such an object has to be considered unlikely. 

Of greater interest are the 0.524 Mq and 0.625 Mq he- 



lium burning models which have three branches in Fig. 3 
(plotted in red and green respectively). Helium burning 
arises because the last thermal pulse occurs in the post- 
AGB phase, transforming the hydrogen burner into a he- 
lium burner; the consequence is rapid evolution back to the 
AGB where post-AGB evolution begins once more as a he- 
lium burner. The longer post-AGB evolution time implied 
by helium burning, through crossing the Toft — logy space 
defined by Fig. 3 three times, makes detection more likely. 
J16007-I-0748 could therefore have a sdO star which is a 
post-AGB helium burning object with M ~ 0.6 Mq; the lu- 
minosity would the be ~ 10^ Lq which was consistent with 
observation as discussed below. 

Radial velocities obtained for the K3 V and sdO stars 
were found to be similar, although with large standard de- 
viations, suggesting they constitute a physical binary. The 
absolute magnitude (Mv) of the K3 V star indicated a dis- 
tance (d) of ~ 3 kpc. Angular radii then gave radii of ~ 0.08 
and ~ 0.7 Rq for the sdO and K3 V stars respectively, with 
a corresponding sdO star luminosity ~ 150 Lq . The Table 1 
logg (sdO) resulted in a sdO star mass of ~ 0.04 Mq. For 
the sdO star to have a mass nearer ~ 0.4 Mq, for reasons 
discussed above, its distance needs to be about a factor of 
three larger. 

A possible associati on of the sdO star with pos t- 
RGB (Driobc ot al. 'l998) or post-AGB ijBlocke j [l995al lbl: 
[Schonbcrncr 1979, 1983) evolutionary tracks would imply 
respective luminosities of ~ 10^'^ and ~ 10^'^ Lq, consis- 
tent with d ~ 9 kpc. A d ~ 9 kpc places J16007-I-0748, 
given its galactic coordinates at ~ 7 kpc from the Galactic 
Centre. For J16007-f 0748 to be a physical binary, the late- 
type companion would then have Mv ~ 4.2 rather than the 
expected Mv 6.7 for a normal K3 V star. 

A higher My for the l ate-type compa n ion co uld be ex- 
plained if the Paper I and IFontaine et al.l ([200^) classifica- 
tions for the late-type companion are correct, and not the 
one presented in this paper; for this to be satisfactory, cor- 
roboration is needed for the implicit Fontaine et al. assump- 
tion of Eb-v ~ 0.0. If the late-type companion had evolved 
off the Main Sequence and is in reality a K3 IV star, the 
higher luminosity and radius might be explained; but this 
is also unsatisfactory because the Section 2 analysis found 
less satisfactory fits to the SDSS energy distribution were 
obtained when lower gravity synthetic spectra were used for 
the late-type companion. Accretion of material from the sdO 
progenitor could result in the K3 V star appearing to be 
over-luminous; a theoretical justification would be needed 
and is beyond the scope of the present paper. 



5 NEW PHOTOMETRY 

New photometry was obtained with the same instrumenta- 
tion and methods used for the observations presented in Pa- 
per I. In brief, observations were made during a "core" run 
of fourteen consecutive nights in 2006 June/ July, followed 
by three consecutive nights one week later and a further two 
nights about a month later. All observations were made with 
the Universit y of C ape Town CCD photometer (UCTCCD; 
IO'Donoghu3ll995l ) on the 1.9-m telescope at the Suther- 
land site of the South African Astronomical Observatory 
(SAAO). The UCTCCD was used in frame-transfer mode 
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Table 2. Observation log for J16007+0748 (1.9-m telescope + 
UCTCCD). 



Date 


JD 


Run 


Comments 


2006 




(hr) 




Jun 28/29 


245 3915 


5.6 




/ ou 




O. I 




ou/ Ul 








Jul 1/2 


3918 


5.5 




2/3 


3919 


5.5 


software crash 


cS/4 


oyzu 


^ fi 
o.u 




4/5 


3921 


1.7 


cirrus 


5/6 


3922 


2.8 


cirrus 


7/8 


3924 


5.1 




8/9 


3925 


3.4 


some cloud 


9/10 


3926 


4.9 


cirrus 


10/11 


3927 


5.2 




11/12 


3928 


5.2 




16/17 


3933 


2.8 




17/18 


3934 


5.0 




18/19 


3935 


4.7 




Aug 17/18 


3965 


2.1 




20/21 


3968 


2.2 





which has the advantage of losing no time on target when 
the CCD reads out, enabUng continuous or "high-speed" ob- 
serving. 

On the 1.9-m telescope, the 22 n pixels of the CCD are 
equivalent to 0.13 arcseconds at the detector, so that it is 
normal to use at least 3x3 prebinning for optimal data ex- 
traction, unless the seeing is better than about 1 arcsecond. 
The first week was almost uniformly good and all observa- 
tions were made with 3x3 pre-binning; the second week 
was less good, though mostly usable, and most of the obser- 
vations were made with 4x4 pre-binning. All integration 
times were 10 seconds, except for the runs on June 28/29 and 
July 4/5, where a time of 12 seconds was used. These values 
were considered a reasonable compromise between obtaining 
as good a signal-to-noise as possible and having adequate 
temporal sampling to resolve rapid variations. Conventional 
procedures (bias subtraction, flat field correction and so on) 
were followed with magnitude extraction being based on the 
po int-spread function of the DoPho t program as described 
bv lSchechter. Mateo fc Sahal ^199^ ). 

All photometry of J16007-I-0748 was corrected differen- 
tially to remove any rapid transparency variations using the 
nearby bright star (about 20 arcseconds SE of J16007+0748; 
see the chart given in fig. 1 of Paper I) . Since, in general, field 
stars will be quite red, we might expect differential extinc- 
tion effects to be significant, and we have further corrected 
for residual atmospheric extinction effect by removing a sec- 
ond order polynomial from each night's observations. This 
means that we might be removing real changes in stellar 
brightness on time scales of a few hours but this has to be 
accepted. 

A log of the "high-speed" or continuous photometry 
obtained for J16007-I-0748 is given in Table 2. (Note that 
runs much longer than about 6 hours are impossible for a 
star at declination ~ -1-8° from Sutherland without reaching 
undesirable air mass values). A sample of a differentially 



Figure 4. 2006 June/July observations for J16007-I-0748. Num- 
bers at the right-hand edge are JD — 2453900. 

corrected light curve is shown in fig. 2 of Paper I; that fight 
curve is very similar to all the photometry analysed here. 

Fig. 4 shows the distribution in time of the 2006 
June/ July observations. It is clear that they fall into two 
main sections, the first six nights and the last five. It can 
also be seen from the figure that the scatter is rather worse 
in the last five nights - a result of somewhat poorer weather 
and seeing conditions, on average. 



6 FREQUENCY ANALYSIS 

The frequency analyses described in this section were carried 
out using software which produces Fourier a mplitude spectra 
following the Fourier transf orm method of iDeemin 3 (|l975l ) 
as modified bv lKurtj l| 19851 ). 

6.1 The "core" data 

The observations shown in Fig. 4 lend themselves in an ob- 
vious way to our preferred practice of reducing the observa- 
tions initially in two approximate "halves". We have taken 
the first six nights and the last five nights and determined 
the Fourier amplitude spectra for these separately. The base- 
lines of the two data sets correspond to resolutions of about 

2.2 and 2.7 /xHz respectively, and the rms scatter in the 
background noise of the Fourier amplitude spectra (avoid- 
ing any obvious signal spikes) are about 0.4 and 0.7 mmag. 
Fig. 5 shows the results after pre-whitening each set by the 
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Table 3. Six frequencies and amplitudes extracted simultane- 
ously plus three frequencies extracted singly (and allowing for 
aliasing) from the analysis represented in Figures 5 and 6. For- 
mal errors in the amplitudes are ~ 0.3 mmag for the first 6 nights 
and 0.6 mmag for the last five nights. 



First 6 nights Last 5 nights 



Freq 


± 


Amp 


Freq 


± 


Amp 


Paper I 






mmag 




/iHz 


mmag 




8379.48 


0.01 


38.1 


8379.53 


0.02 


38.0 


/i 


9091.32 


0.05 


5.2 


9091.08 


0.11 


6.5 


/2 


16758.91 


0.06 


4.2 


16759.00 


0.15 


4.7 


fa 


14186.28 


0.05 


4.9 


14186.86 


0.16 


4.3 


U 


13073.82 


0.08 


3.4 


13074.67 


0.10 


7.0 


/7//13 


8461.7 


0.12 


2.2 


8462.4 


0.25 


2.8 


/5 


9663.0 


0.14 


2.1 


9686.2 


0.34 


2.4 


/6 ? 


8506.3 


0.14 


1.9 


8510.5 


0.33 


2.1 


hi 


15939.3 


0.23 


1.3 


15940.4 


0.35 


2.3 


fe 



rnm&g 



6 


1 




1 

(a) _ 


1 
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Figure 5. Fourier amplitude spectra for (a) the first six nights 
and (b) the last five nights (see Fig. 4). Both have been pre- 
whitened by the very strong frequency near 8380 fiHz, better to 
display the weaker frequencies. The corresponding spectral win- 
dows are plotted at the right with the same frequency scale and 
normalised amplitudes. 

rmnag 




8000 12000 16000 20000 



Figure 6. Fourier amplitude spectra for (a) the first six nights 
and (b) the last five nights (see Fig. 4) after the removal of six 
common frequencies. 

6.2 The first six "core" nights 

Because amplitude variation - whatever the cause - will 
confuse the analysis, we have taken the first six nights (which 
are clearly of better quality than the rest of the observations) 
and carried out a similar process to that in Section 6.1, in 
this case splitting the observations into two sets of three 
nights. The resolution for each set is about 5 /iHz and the 



dominant frequency near 8380 //Hz to make the weaker fre- 
quencies more visible. The difference in background noise is 
also clear in that figure and is a result of the poorer quality 
of the second half of the observations (See Fig. 4). 

We have initially extracted frequencies from the ampli- 
tude spectra one at a time and then - for several frequen- 
cies common to both sets of obs ervations - have carried out 
simultaneous least-squares fits (|Deemind IToei ) to the ex- 
tracted frequencies, followed by a search for any remaining 
common frequencies. By continuing this procedure we find 
six well-established frequencies which are listed in Table 3. 

The Fourier amplitude spectra of the observations with 
six frequencies removed are shown in Fig. 6, in which it 
appears that there might still remain common frequencies. 
Such were not immediately obvious in our initial search but, 
allowing for 1 or 2 cycle/day aliasing (corresponding to dif- 
ferences of 11.6 and 23.1 /iHz), we can find three more com- 
mon frequencies which are also listed in Table 3. The final 
column of that table lists the frequency designations, fi, 
given in Paper 1. 

Using Paper I nomenclature, we easily recover /i to /s; 
fe - 9651.9 fiHz in Paper 1 - is almost certainly the peak 
we find at 9663.0 /iHz with a one cycle d"'^ alias either in 
Paper I or here. The peak near 13074 /iHz is probably fr 
with a two cycle d~^ alias, but could also be /13 (13072 /xHz 
in Paper I). Somewhat surprisingly, /g is a marginal detec- 
tion in the current work and we do not recover /g - /12 
even at only twice the background noise in the present data, 
whereas they were detected in Paper I at four to six times 
the background. Curiously, /14 does seem to be present and 
at about the same amplitude as in Paper I. 

In a sense, this is a disappointing result - with far more 
observations over a slightly longer baseline than the Paper I 
data, we had expected to recover all the original frequencies 
and hoped to find even more. The current result probably 
reflects the fact that some of the frequencies have variable 
amplitude; this was evident already in Paper 1 (see fig. 3 
of that paper) and is also suggested in our Fig. 5 and Ta- 
ble 3 where the peak near 13074 ^Hz doubles in amplitude 
between the first six nights and the last five; in a similar 
comparison, the other peaks have constant amplitude within 
about 2 (J to 3 (J. 
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background noise is around 0.5 mmag. This subset of the 
observations is comparable with that of Paper I - almost 
30% more data but obtained over a baseline of six nights 
rather than ten. 

After the removal of eight common frequencies with am- 
plitudes greater than about four times the background, there 
still seem to be significant peaks in the amplitude spectra 
(Fig. 7). If, as before, we allow for aliasing and also search 
down to three times the background (given that we have 
two separate sets of data), five more common frequencies 
are found. Again, we recover /i to and /14, all at four 
times the background or better (see Table 4). In addition, 
we have weak evidence for /g, /lo and /12 and maybe two 
new frequencies, but these results are not totally convinc- 
ing. We still find no significant evidence for the presence of 
/g and /ii; this does not mean we doubt their existence in 
the Paper I data, indeed, they seem well-established there. 
Rather, it seems more likely that we are seeing more evi- 
dence for amplitude variation. 



6.3 The June/ July and August data 

Analysing separately the three nights from mid- July (16/17 
- 19/20), we easily recover the frequencies /i to /t, well 
above the background noise (allowing for some one- and 
two-cycle/day aliasing). However, when we add these ob- 
servations to the core data, the analysis is less clear. The 
strongest frequencies are still recovered but appear to be 
accompanied by weak peaks within about a ^Hz or so. This 
immediately suggests that variable amplitudes might be the 
cause (whether due to unresolved pairs of frequencies or ac- 
tual amplitude variation in a single mode) and that these 
are made clearer by the longer baseline. But it is also possi- 
ble, given that J 16007-1-0748 is a spectroscopic binary of un- 
known orbital period (see Paper I), that the orbital motion 
is creating an effective phase shift in the later results rela- 
tive to the earlier. In this case, any derived frequency would 
be affected and extraction of a sine wave with a slightly in- 
correct frequency could result in cyclic residuals in the data 
which would then be extracted as a nearby, much weaker 
signal. Detailed spectroscopic determination of the orbital 
period would help to resolve this issue - and might even 
enable the pulsation phase shifts to be used to measure the 
size of the projected orbit by determining th e light trave l 
time. This was attempted by iKilkennv et all l| 19981 , 120031 ) 
for the eclipsing, pulsating PG 1336-018, but resulted only 
in upper limits of about a light second for the projected or- 
bital diameter - but that star has a period of only about 
0.1 d and the problem would be more tractable for a longer 
period binary. 

The August observations were analysed separately and 
revealed frequencies /i to /s but both series of observations 
were short (around 2 hours) and are separated by a month 
from the other observations, so we have not attempted an 
analysis of the merged data. 
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Figure 7. Fourier amplitude spectra for (a) the first three nights 
and (b) the second three nights (see Fig. 4) after the removal of 
eight common frequencies. 



Table 4. Eight frequencies and ampHtudes extracted simulta- 
neously from the first six nights split into two halves, plus five 
frequencies extracted singly (and allowing for aliasing). Formal 
errors in the amplitudes are about 0.3 mmag for both halves. 



First 3 nights 



Second 3 nights 



Freq 


± 


Amp 


Freq 


± 


Amp 


Paper I 


/iHz 




mmag 


/iHz 


/xHz 


mmag 




8379.48 


0.02 


38.3 


8379.50 


0.02 


38.2 


/i 


9091.8 


0.18 


4.1 


9090.8 


0.12 


6.4 


/2 


16759.0 


0.16 


4.7 


16758.9 


0.18 


4.1 


h 


14186.3 


0.15 


5.0 


14186.1 


0.16 


4.7 


h 


13074.2 


0.16 


4.6 


13072.4 


0.29 


2.6 


h 


8461.1 


0.42 


1.8 


8460.4 


0.24 


3.2 


h 


9664.1 


0.26 


2.9 


9662.9 


0.40 


1.9 


k 


8518.1 


0.40 


1.9 


8516.8 


0.36 


2.1 


hi 


12816.5 


0.39 


2.4 


12815.2 


0.70 


1.3 




12521.3 


0.37 


2.0 


12536.0 


0.33 


2.2 


fl2 ? 


17483.9 


0.51 


1.5 


17459.0 


0.38 


1.9 




13667.6 


0.50 


1.5 


13659.6 


0.41 


1.8 


flO ? 


15921.6 


0.56 


1.7 


15939.7 


0.71 


1.3 


/8 



6.4 Summary 

The frequency analysis of the observations reported here has 
found the same frequencies /i to fj as found in Paper I. It 
is possible, we have also detected /g and /14 and maybe 



others, though these detections are much less secure. It ap- 
pears, somewhat paradoxically, that shorter data sets yield 
"better" results and this is likely to be due to amplitude 
variation (or equivalently, frequency beating) or the effects 
of the binary orbital motion - and probably both of those. 
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Figure 8. Evolutionary tracks and predicted positions of 
J16007+0748 in a T(,ff-log g diagram. Diamonds mark sdO models 
along the tracks analysed in search for instability. The analysis of 
those marked with a filled diamond is explicitly shown in the pa- 
per (see text for further details). The large upper rectangle shows 
the log g-Tefl and its uncertainty as determined in the present pa- 
per. T he lower and smaller r ectangle is the corresponding position 
which iFontaine et al. determine. 



7 PULSATION MODELLING 

A pulsation analy sis of structural models of 

iRodrfeuez-Lopez et all (|2008l fl was carried out within 
the logg-Teff 68 per cent confidence interval determination 
for J16007-h0748. 

To create the starting points for sequences labeled 1 
through 5, a star of initial mass 1 Mq and initial heavy 
element abundance Z = 0.02 was evolved from the pre- 
Main Sequence stage to the Red Giant Branch (RGB), 
through the helium core flash and to the Zero Age Horizon- 
tal BranchJ^AlIB)^Tlie cod^^ JMSTAR, is described 
bv lLawlor fc MacDonaldl dioosl ) and references therein. In 
order to produce models of hot subdwarfs l|D'Cruz et al.l 
1 19961 ). the mass-loss rate on the RGB is enhanced compared 
to typical observed values, and it is characterised by a value 
of the Reimer's mass-loss parameter riR = 0.675. For com- 
parison, typical observed RGB mass-loss rates correspond 
to riR = 0.35 - 0.4. Because of the high mass-loss rates, the 
star leaves the Giant Branch before the helium core flash oc- 
curs. The mass of the hydrogen- rich envelope, ~ 0.01 A/0, at 
this point is sufBciently small so that the star arrives on the 
ZAHB relatively hot, T^ft ~ 20 000 K, and has the charac- 
teristics of a subdwarf-B star. As a greater than solar heavy 
element abundance in the driving zone seems necessary for 
pulsational instability, an ad hoc extra mixing is introduced 
for the evolution beyond the ZAHB. In the stellar evolution 
code, convective mixing is modelled by use of diffusion equa- 
tions, with the turbulent diffusivity calculated from mixing 
length theory. Other mixing processes can be included by 



^ The evolutionary s e quenc es in this paper correspond to 
iRodn'guez-Lopez et al.l 1I2OO8I ) sequences 10, 11, 12, 13, 14 and 
7, respectively. Models 1.1, 3, 4 and 5 in the table correspond, 
respectively, to Models 10.3, 12, 13 and 14 of the same paper. 



Figure 9. Mass fractions of the main chemical constituents of 
representative Models; model 1 (blue), 2.1 (green), 5.1 (orange) 
and 6 (magenta). X(Mg) refers to the mass fractions of all the 
other elements. Top right numbers designate panels' numbers (see 
text). 
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Figure 10. Main physical parameters of representative models 
(colour key the same as in the other figures). Note that Panels 1 
and 2 show only the outer envelope, while Panels 4 and 5 show 
only the innermost regions of the star. 



adding terms to the turbulent diffusivity. For our ad hoc 
mixing, the turbulent diffusivity is 

1 rjMT 



^ad hoc — 



(1) 



(47rr2p)^ ^ 

where rj is the mass in a zone and Mt the total mass 
of the star. The degree of extra mixing is determined by 
the ratio of the input time scale r to the evolution time 
scale. A consequence of extra mixing is that core material is 
transported to the surface, increasing heavy element abun- 
dances in the hydrogen-rich envelope. Also, due to mixing of 
helium to higher temperatures, the star experiences helium 
shell flashes that drive convective mixing of heavy elements 
to the surface further increasing Z. The starting points for 
Sequences 1 and 2 are taken from the post-ZAHB track when 
Tcfi = 45 000 K. Similarly the starting points for Sequences 
3, 4 and 5 are taken from the post-ZAHB track when T^s = 
70 000 K. Extra mixing is turned off for the evolution from 
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Table 5. Effective temperature, surface gravity and envelope mass fractions of the equilibrium models. The column labeled X(Mg) refers 
to the mass fractions of all the other elements. Z is the current metallicity of the model. The capital letter beside some model's number 
refers to their plotting colours. 



Model 




log g 


X{H) 


X{He3) 


X(He4) 


X{C) 


X{N) 


X(0) 


X(Mg) 


Z 




(K) 




















1 B 


65 900 


5.04 


0.20 


2.5E-05 


0.64 


7.0E-02 


l.OE-02 


6.2E-02 


1.9E-02 


0.16 


1.1 


69 000 


5.13 


0.19 


2.4E-05 


0.65 


7.1E-02 


l.lE-02 


6.2E-02 


1.9E-02 


0.16 


1.2 


72 400 


5.25 


0.20 


2.2E-05 


0.65 


7.2E-02 


l.lE-02 


6.3E-02 


1.9E-02 


0.15 


1.3 


74 700 


5.33 


0.16 


1.5E-05 


0.67 


7.6E-02 


l.lE-02 


6.7E-02 


2 nF,-n2 


0.17 


2 


66 500 


5.13 


0.31 


8.4E-05 


0.57 


5.2E-02 


8.7E-03 


4.8E-02 


1.5E-02 


0.12 


9 1 

Z. 1 


70 son 

/ u OUU 


o.ou 


n 9Q 


/ .orlj-UO 


U.Oo 


r: 09 

o.orij-uz 


o.yrlj-uo 




1 . 0-Cj-UZ 


u. 10 


2.2 


72 700 


5.37 


0.26 


5.3E-05 


0.60 


6.0E-02 


9.4E-03 


5.4E-02 


1.7E-02 


0.14 


2.3 


75 000 


5.43 


0.26 


5.0E-05 


0.60 


6.1E-02 


9.5E-03 


5.4E-02 


1.7E-02 


0.14 


3 


70 000 


5.16 


0.19 


2.3E-05 


0.65 


7.1E-02 


l.lE-02 


6.3E-02 


1.9E-02 


0.16 


3.1 


71900 


5.22 


0.19 


2.2E-05 


0.65 


7.2E-02 


l.lE-02 


6.3E-02 


1.9E-02 


0.16 


3.2 


73 600 


5.30 


0.16 


1.5E-05 


0.67 


7.6E-02 


l.lE-02 


6.6E-02 


2.0E-02 


0.17 


3.3 


75 000 


5.34 


0.16 


1.5E-05 


0.67 


7.6E-02 


l.lE-02 


6.7E-02 


2.0E-02 


0.17 


4 


70 000 


5.20 


0.16 


1.6E-05 


0.66 


7.6E-02 


l.lE-02 


6.6E-02 


2.0E-02 


0.18 


4.1 


74 400 


5.32 


0.16 


1.5E-05 


0.67 


7.6E-02 


l.lE-02 


6.7E-02 


2.0E-02 


0.18 


5 


70 000 


5.43 


0.48 


3.0E-04 


0.45 


2.8E-02 


6.0E-03 


2.9E-02 


l.OE-02 


0.07 


5.1 O 


73 600 


5.54 


0.43 


2.2E-04 


0.48 


3.4E-02 


6.7E-03 


3.4E-02 


l.lE-02 


0.09 


6 M 


70 800 


5.98 


0.59 


1.9E-04 


0.36 


7.3E-03 


4.3E-03 


2.4E-02 


1.5E-02 


0.05 



these starting points. For Sequences 2 and 5 the mass-loss 
rate is set to zero. For Sequences 1, 3 and 4 mass- los s from a 
line-d riven wind is assumed with the rate taken from lAbbottI 
l| 19821 ) but scaled in proportion to Z"'^. All models at this 
point of evolution have mass 0.479 Mq. 

The starting point for Sequence 6 is from the track for 
a star of initial mass 1 Mq, initial metallicity Z = 0.05, 
and riR = 0.6. Since the heavy element abundance is greater 
than solar, no extra mixing was used. The evolution of this 
model star is more normal in that it is the helium core flash 
that terminates the RGB. At the ZAHB, the hydrogen enve- 
lope has a mass ~ 0.06 Mq and hence the star is cool with 
Teff ~ 4300 K. However, due to mass-loss primarily during 
the helium shell burning phase, the mass of the hydrogen en- 
velope is sufficiently reduced to avoid the AGB phase. The 
starting point from Sequence 6 is taken at Tctt — 45 000 K 
as the star is evolving at roughly constant luminosity to the 
blue, on its way to becoming a white dwarf. At this point 
the stellar mass is 0.497 Mq. 

The evolutionary tracks described by each evolutionary 
sequence in the log g-T^s diagram are shown in Fig. [H] We 
note that Sequences 1, 3 and 4, describe the same evolution- 
ary track, although equivalent models have different enve- 
lope compositions. A stability analysis was carried out for 
models marked with diamonds. Below, the analysis is shown 
for representative Models 1, 2.1 and 5.1 (marked with filled 
diamonds) located around the centre, and two of the corners 
of the box. In additio n, Sequence 6 which c rosses the log g- 
TbS determination of iFontaine et al"] (I2OO8I ) is also plotted 
in Fig. [S] and a model from this sequence within their spec- 
troscopic box analysed for comparison. 

Table Oshows the models' Tofi, \ogg and envelope mass 
fractions, which, for the selected models (marked with the 



initial of their plotting colour) , are plotted in Fig. |^. Panel 
1 shows that the cores of all models are devoid of hydrogen 
inside \ogq ~ —2.5, where hydrogen shell burning begins. 
The models are also devoid of helium in the core, until it 
starts to burn in a shell around logg ~ — 1. Panel 2 shows 
He'' which is an intermediate product of the proton-proton 
chain reaction (PP I) produced in the hydrogen burning 
which creates He* nuclei. He^ can also directly react with 
He* nuclei and, via intermediate processes, produce two He* 
nuclei (PP H process). Panels 3 and 5 shows that the core 
of the models is composed of the residuals of the helium 
burning: carbon, produced in the 3-a process, and oxygen, 
as a result of the carbon capturing helium nuclei. N^*, in 
Panel 4, is produced as a subproduct of the CNO cycle (C^^ 
recombining with protons, shown in the decrease of carbon 
at logg ~ —2) which creates more C'^^ and He*. The last 
panel display the mass fraction of all the remaining heavy 
elements. 

Some representative physical parameters of the selected 
models are shown in Fig. 1101 Panels 1 and 2 display the ex- 
pected behaviour of pressure and density, respectively. Panel 
3 shows the logarithm of the opacity: the largest opacity 
bump corresponds to the heavy elements partial ionisation 
zone, and the secondary bump at logg ~ —6 to the C/0 
partial ionisation zone. We note that for Model 6 (magenta) 
the opacity is larger, although the derivative of the opacity - 
implied in the k mechanism driving process- is qualitatively 
similar. Panel 4 shows the nuclear energy generation rate 
normalised to its maximum: the deepest peak corresponds 
to the helium burning shell, while the shallowest is that of 
the hydrogen burning shell. We note that for Model 6 the 



^ All figures are plotted as a function of the fractional mass depth 
log q = log(l — Mr/Mj-), which better samples the outer envelope 
of the star. 
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Figure 11. LeftiGrowth rates of representative models. Asterisks mark modes g76, gl5 and p2 of Model 5.1, representative of its L, I 
and H regions, respectively. The two vertical dashes delimit the range of instability found for J 16007+0748. Right: Growth rates of two 
models with same log g and different Tgff . (See text for details) . 
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Figure 12. Left, Center & Right: Energy and opacity for the quadrupole modes of Model 5.1 g76, 51I5 and p2, belonging to L, I and H 
zones, respectively (see text for details). All plots have been scaled to arbitrary units. 



contribution of the helium burning is very low compared 
to that of the hydrogen; and that hydrogen burning takes 
place deeper in the rest of the models. Panel 5 shows the 
luminosity normalised at the surface of the star with the 
contributions of both burning shells; for all the models, ex- 
cept Model 6, the luminosity achieves negative values at the 
location of the He burning shell. This is due to the increase 
in the helium burning rate produced by the extra-mixing of 
He, C, N and O in the envelope. The energy so produced 
diffuses inwards to the semi-degenerate material right be- 
low the burning shell, that expands and cools producing the 
negative luminosities. 

Finally, Panel 6 shows on one hand, the logarithm of 
Brunt- Vaisala frequency (solid lines) which shows local max- 
ima at the zones with maximum composition gradients. The 
deepest and largest maximum is due to the transition be- 
tween the C/0 core and the He burning shell, while the 
secondary peak -stronger for Model 6, as the transition is 
more abrupt (see Panel 1 in Fig. |9])- is due to the He/H 
transition. All the energy transport in the star is done by 
radiation except at the location of the heavy elements par- 
tial ionisation zone where A^'^ < 0, meaning that the energy 
transport is more efficiently done by convection. Panel 6 also 
shows the Lamb frequency (dashed lines) for I = 2 modes. 

The models so constructed were subject to a non adi- 
abatic stability analysis carried out with GraCo non adia- 



batic code of oscillations (|Mova fc Garridd |2008| ). Modes 
were calculated from / = to i = 4 in the 600- 
20 OOP I-lHz frequency range, and the normalized growth rate 
rj l|Dziembowski et al.lll993r ). was used as an indicator of 
stability. In our convention, positive (negative) values of the 
growth rate mean the mode is unstable (stable). Fig. llll (left') 
shows the growth rate values of each of the selected mod- 
els, all of which were found stable in the analysed frequency 
range. 

Each model presents three different regions according 
to the behaviour of the growth rate. Each of these regions 
cannot be delimited in frequency, as they are model depen- 
dent, but it is useful to label them generally as low (L), 
intermediate (I) and high frequency region (H) to describe 
their general behaviour: 

• The very low frequency region (L), with a slight ten- 
dency to instability, comprises more modes (extending to 
higher frequencies) and its driving is increased as the model 
logg increases, and it is only slightly temperature depen- 
dent. 

• The intermediate frequency region (I) is defined by 
growth rate values of 77 = —1, indicating its high stability. 
This region also shifts to higher frequencies and includes 
more modes as the log g of the models increases. 

• The high frequency region (H) corresponds to frequen- 
cies over the I region which have a certain tendency to in- 
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stability. This region also shifts to higher frequencies and 
comprises a wider range in frequencies as the model log 17 
increases. Besides, when the model T^s increases at constant 
log g, there is an overall increase in the tendency towards in- 
stability (see Fig. [11] right). We also verified that, for models 
with equal metallicity, those with higher T^s had somewhat 
higher growth rate values. 

The described behaviour of the L, I and H regions is also 
found for models within the same evolutionary sequence. Be- 
sides, as expected, the frequency of the fundamental mode 
increases when the log 3, or mean density of the model in- 
creases. Therefore, the high frequency region is associated 
with low-radial order p-modes for the models with lower 
logp, and shifts into the p-mode region as the model logg 
increases. 

We found Model 5.1 to be most favoured for instabil- 
ity within the range of pulsation frequencies observed for 
J16007-I-0748 (roughly 8 to 16 mHz), this has the highest 
log g and opacity among our evolution sequences. 

We present a stability analysis carried out for Model 
5.1 as representative of all the other extra-mixed models. 
We examine the local energy exchange for modes belong- 
ing to L, I and H zones, respectively, through the derivative 
of the work integral, dW / dlogq, whose negative (positive) 
values indicate local contribution to driving (damping) of 
the mode. Fig. 1121 shows the energy balance for representa- 
tive modes glQ, 17I5 and p2, with the opacity plotted in the 
right axis for a better understanding of the plot. Mode glQ 
(left) is representative of the low frequency range, compris- 
ing from 0.6 to ~2.5 mHz, where modes achieve a certain 
destabilisation. This is caused by the existence of a driving 
region at the location of the C/0 partial ionisation zone at 
log q ~ —6. Mode gib (centre) belongs to the intermediate 
frequency range, with cyclic frequencies between about 2.5 
to 6 mHz, where modes are highly stable. This is explained 
as at these frequencies, the maximum energy interchange 
occurs in a damping region. Mode p2 (right) is representa- 
tive of the high frequency range, from 6 mHz and above. 
The growth rate increases, and so the tendency to driving, 
due to the existence of a driving region at the location of 
the Z-bump at logg ~ —10.5 for this model. The increase or 
decrease of the driving and damping regions with frequency, 
around the Z-bump, accounts for the oscillating profile of 
the growth rate in this region. 

By comparison. Model 6, belonging to Sequenc e 6 which 
passes through the Xogg-T^s region which Fonta irie et al.l 
l|2008l ) obtain, presents only a destabilisation for very low 
and very high frequencies (in the range 20 to 50 mHz, out of 
the range of the plot), due to the C/O and heavy elements 
partial ionisation zones, respectively. The observed instabil- 
ity frequency region for J 16007-1-0748 is highly stable for 
this model, as any significant energy dampes the modes. 



7.1 Discussion 

Our pulsational study has to be compared with the dis- 
covery of the mechan ism of oscillations of J 16007-1-0748 by 
iFontaine et al.l l|2008l ). These authors achieve real driving 
of modes within the instability range observed, using non- 
uniform metallicity models as a result of the diffusive equilib- 
rium between radiative levitation and gravitational settling. 



The need for a non-uniform distribution of iron with depth 
to achieve an accurate description of the instability strip 
for sdBs had already been established by the same group 
dCharpinet et aQ 1997h . a nd it has now been confirmed by 
iFontaine et al. I 2008h and ICharpinet. Fontaine &: BrassardI 
l|2009l ) to be a crucial ingredient also to achieve real driv- 
ing in sdO models. All of our models have been built with 
uniform metallicity, and hence our inability to achieve real 
driving of the modes. 

However, from several tests carried out in our pulsation 
analysis we still can conclude that for sdO models within 
our derived spectroscopic box, those with log g values higher 
than 5.30 are more likely to be driven in the observed fre- 
quency range of J16007-I-0748. Besides, for models with the 
same logg, hotter models have a higher tendency to drive 
pulsations. Hence, the pulsation analysis favours the higher 
effective temperature and higher logp estimate within our 
spectroscopic box. Among the analysed models. Model 5.1 
with Tcff=73 600 K and logp = 5.54 was found the most 
favoured for driving within the unstable frequency range ob- 
served for J16007-f0748. 



8 CONCLUSIONS 

We have presented a spectroscopic analysis of J16007-I-0748, 
a binary system containing the only pulsating sdO star 
known to date, and a late type companion. A flux calibrated 
spectrum from the sixth data releas e of the SPS S was used 
for the analysis. We used tlusty (jHubenv fc L anz 199^) 
to compute a grid of non-LTE model stellar atmospheres 
for the sdO co mprising H a nd He, which were further used 
with SYNSPEC (IZborillll996l . and references therein) to com- 
pute a grid of non-LTE syntheti c spectra. For t he lato-type 
companion, synthetic spectra bv lMartins et al] ((2005) were 
used, and a null projected rotational velocity was assumed 
for both stars. A genetic algorithm optimisation was used 
to estimate parameters for both binary components given in 
Table 1. 

We co mpared our results w ith atmospheric parameters 
derived bv IFontaine et al] (|2008l ) and find that our Toff esti- 
mate is consistent, although with larger error limits; and our 
log(? is lower than theirs. As for the cool companion type, 
the difference in spectral type estimate may arise from the 
fact that we have included an interstellar reddening factor, 
which required a redder model for the companion. 

A discussion of possible evolution scenarios for the 
sdO component took our results into account. The de- 
rived atmospheric parameters classify J160074-0748 as a 
helium-rich sdO, although showing Tcff and \ogg values 
above those canonic ally considered for these objects (see i.e. 
IStroeer et ai]|2007l ). We considered J16007-I-0748 being a 
post-RGB object fo llowing an evolutionary track of 0.4 Mq 
l|Driebe et al.|[l998l ), but it is not clear how the helium en- 
richment would have occurred. We also cons i dered t he ob- 
ject to be a post-AGB st ar following iBlockeij (|l995al lbh and 
^Schonberned (|l979l . Il983l ) evolutionary tracks. For most of 
the tracks the rapid evolution from the AGB to the pre- 
white dwarf stage renders the observational detection un- 
likely. However, two 'born-again' AGB tracks from Blocker, 
with 0.5 and 0.6 Mq, having longer evolutionary time, make 
the observational detection more likely. The possibility of the 



The sdO pulsator J16007+0748 13 



system to be physically bounded, as it seems to be implied 
in the derived radial velocities, is also discussed. However, 
uncertainties in the derived absolute magnitudes and stellar 
masses prevent us from a final conclusion. 

We presented new and more extensive photometry gath- 
ered during nineteen nights on the 1.9-m telescope at SAAO. 
The frequency analysis recovers with certainty seven out of 
ten frequencies established in an earlier paper. We conclude 
that the shorter data sets give better results, possibly due 
to amplitude variations that may be caused by the binary 
orbital motion, which needs further investigation. Therefore, 
a spectroscopic determination of the orbital period is very 
desirable. 

Finally, a pulsation analysis of uniform-metallicity sdO 
models within the atmospheric parameter range showed no 
modes were actually excited, but useful conclusions were 
obtained. We found models with \ogg > 5.30 were more 
likely to be unstable at the frequency range observed for 
J16007-I-0748. Besides, for models with the same logg, 
higher effective temperatures increase the overall driving 
rate. Therefore, among our models, the most favoured for 
instability within the observed frequency range has atmo- 
spheric pa rameters Tes = 73 600 and log g = 5.54. We re- 
mark that iFont aine et al] l|2008h achieve actual driving of 
p-modes within the instability range of J 16007-1-0748 with 
models including radiative levitation, which they identify as 
key ingredient in driving the oscillations. 
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